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Abstract 

Background: Most of the natural variation in flowering time in Arabidopsis thaliana can be attributed to allelic 
variation at the gene FRIGIDA {FRI, AT4G00650), which activates expression of the floral repressor FLOWERING LOCUS 
C [FLC, AT5G10140). Usually, late-flowering accessions carry functional FRI alleles [FRI-wt], whereas early flowering 
accessions contain non-functional alleles. The two most frequent alleles found in early flowering accessions are the 
ones present in the commonly used lab strains Columbia (FRI-Col) and Landsberg erecta (FRI-Lef), which contain a 
premature stop codon and a deletion of the start codon respectively. 

Results: Analysis of flowering time data from various Arabidopsis natural accessions indicated that the FRI-Ler allele 
retains some functionality. We generated transgenic lines carrying the FRI-Col or FRI-Ler allele in order to compare 
their effect on flowering time, vernalization response and FLC expression in the same genetic background. We 
characterize their modes of regulation through allele-specific expression and their relevance in nature through 
re-analysis of published datasets. We demonstrate that the FRI-Ler allele induces FLC expression, delays flowering 
time and confers sensitivity to vernalization in contrast to the true null FRI-Col allele. Nevertheless, the FRI-Ler allele 
revealed a weaker effect when compared to the fully functional FRI-wt allele, mainly due to reduced expression. 

Conclusions: The present study defines for the first time the existence of a new class of Arabidopsis accessions 
with an intermediate phenotype between slow and rapid cycling types. Although using available data from a 
common garden experiment we cannot observe fitness differences between accessions carrying the FRI-Ler or the 
FRI-Col allele, the phenotypic changes observed in the lab suggest that variation in these alleles could play a role in 
adaptation to specific natural environments. 

Keywords: Arabidopsis thaliana, Flowering time, FRIGIDA, Vernalization, Natural variation. Allelic series 



Background 

As plants are sessile organisms, adaptation to the environ- 
ment is essential for their survival and reproductive success. 
Mechanisms regulating the response to environmental cues 
enable a proper timing of key events in a plant s life. One 
crucial event, resulting from the integration of endogenous 
and environmental signals, is the switch from vegetative to 
reproductive development. The annual species Arabidopsis 
thaliana occurs in the northern hemisphere in a broad 
range of latitudes differing substantially in day length, 
temperature and other ecological factors [1]. As a result of 
adaptation to specific habitats, Arabidopsis accessions have 
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evolved two main life history strategies. Winter-annuals 
germinate in autumn, survive winter as a rosette and flower 
in the following summer, whereas summer- annuals germin- 
ate in spring or summer and finish their reproduction cycle 
in the same year. Variation between these distinct strategies 
has been associated with allelic variation at the genes 
FRIGIDA {FRI) and FLOWERING LOCUS C (FLC), which 
act epistatically to regulate flowering time [2,3]. Arabidop- 
sis individuals containing functional alleles at these two 
loci flower very late or not at all, unless they receive a pro- 
longed exposure to cold (vernalization). 

FLC encodes a MADS -box transcription factor that 
binds to the promoters of floral initiators such as FLOW- 
ERING LOCUS T (FT) and SUPPRESSOR OF CONSTANS 
OVEREXPRESSION I (SOCI) to repress their transcrip- 
tion. Flowering occurs when FLC is downregulated by pro- 
teins of the vernalization and/or autonomous pathways. 
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reviewed in [4-7]. Transcriptional activation of FLC re- 
quires integrated activity of diverse chromatin remodeling 
and histone-modifying complexes, reviewed in [8]. When 
Arabidopsis plants are vernaUzed, expression of FLC is 
decreased and maintained at reduced level by different 
epigenetic marks in a Polycomb-mediated process in- 
volving long non-coding RNAs [9-11]. The epigenetic 
silencing of FLC is quantitatively modulated and under- 
lies Arabidopsis natural variation for vernalization re- 
sponse [12-14]. Two main haplogroups of FLC have 
been defined mainly by polymorphisms within the first 
intron, an important region to maintain silencing in- 
duced by vernalization [15,16]. These haplogroups were 
shown to underlie differences in flowering time among 
natural accessions of Arabidopsis, but only when a 
functional FRI allele is present [15]. 

Despite the central role of FLC, most of the variation 
in flowering time has been found to correlate with natural 
allelic diversity at FRI [3]. FRI is the founding member 
of a family of seven Arabidopsis proteins that contain 
two coiled-coil domains and show no homology to other 
proteins [2,17]. In order to regulate flowering time, FRI 
builds the scaffold protein of a transcription activator 
complex that mediates diverse chromatin modifications at 
FLC [18,19]. Furthermore, FRI is suggested to be involved 
in co-transcriptional processes that link the function of 
5 ' end capping with transcription and efficient splicing 
of FLC [20]. Vernalization abolishes the effect of FRI 
and silences FLC as described above [9,12]. 

A considerable number of different FRI haplotypes 
have been identified within accessions from a wide range 
of latitudes [14,21,22] or more restricted geographic 
regions [23,24]. Studies on these alleles led to the con- 
clusion that early flowering types evolved by multiple 
independent mutational events from winter-annuals 
containing an ancestral functional FRI allele (FRI-wt) 
[25]. Two distinct deletions in FRI are believed to confer 
early flowering in most of the rapid cycling accessions. 
The Columbia allele {FRI-Col) carries a 16 bp-deletion 
resulting in a premature stop codon and, thus, a truncated 
protein missing a part of the C-terminal [2,14]. But the 
most frequent deleterious FRI mutation in nature is a 
376 bp-deletion combined with a 31 bp-insertion in the 
promoter as observed in Landsberg erecta {FRI-Ler) 
[14,21,22]. This mutation disrupts the translational start 
but, due to a second alternative start codon, a short 
out-of- frame protein might be built [2]. The loss-of- 
function FRI alleles found in Ler and Col are widely 
used as examples of positive selection towards rapid 
cycling accessions [26,27]. 

Although the FRI-Ler and FRI-Col alleles are always 
classified as a single non-functional group, there are evi- 
dences for differences in their effects. First, accessions 
carrying the Ler-type deletion but not the Col-type 



deletion show considerable variation in flowering time 
and FLC levels [14,28]. Then, several studies find vari- 
ation in flowering time associated with the chromosomal 
region of FRI in mapping populations derived from ac- 
cessions containing the FRI-Ler allele crossed to acces- 
sions with a loss-of- function allele [28-30]. In all these 
cases, the Ler allele was associated with a delay in flow- 
ering time. This delay was attributed to additional loci in 
the same region as FRI, although these have never been 
identified. 

This article provides the first robust evidence that the 
Ler allele of FRI is functional. In contrast to the true null 
FRI-Col allele, the FRI-Ler allele is able to induce FLC 
expression, resulting in delayed flowering and increased 
vernalization sensitivity. Nevertheless, the FRI-Ler allele 
has a weaker effect than the fully functional FRI-wt 
allele. The reduced functionality of the FRI-Ler allele 
mainly results from its lower expression putatively due 
to cis regulatory polymorphisms in its promoter. Our 
finding defines a new functional class of FRI alleles, in 
which accessions carrying the FRI-Ler allele would flower 
in between the late and early flowering groups defined 
so far. 

Results 

Effects of the FRI-Ler and FRI-Col alleles on flowering time 
and vernalization response 

We investigated if the FRI alleles present in Col-0 and 
Ler show evidences of different functionality by analyzing 
flowering time data from a recombinant inbred line (RIL) 
population derived from a cross between the two acces- 
sions [31]. Because the effect of FRI on flowering time de- 
pends on the alleles at FLC, we took into account the 
genotype of the RILs at both loci [15]. While Col-0 har- 
bors a functional FLC allele, Ler has been described to 
contain a weak allele with a transposon-like insertion in 
the first intron [25,32,33]. As shown in Figure la, plants 
carrying the FRI-Ler allele flower significantly later than 
plants carrying the FRI-Col allele, although only in the 
presence of the strong FLC allele of Col-0 (two-way 
ANOVA, FRI p = 0.0095, FLC p = 0.0138, interaction 
p = 0.0824). This epistatic interaction suggests the exist- 
ence of a functional FRI allele in Ler, as no other gene 
located in this region of chromosome 4 has been shown to 
delay flowering time through interactions with a locus in 
the region of chromosome 5 containing FLC 

Because strong vernalization requirements are associ- 
ated with functional FRI alleles, we decided to analyze 
variation in vernalization sensitivity among Arabidopsis 
accessions [34,35]. We calculated the decrease in flower- 
ing time in response to vernalization in a published data- 
set containing 126 Arabidopsis accessions classified as 
FRI-Col, FRI-Ler or FRI-wt [22]. Accessions carrying the 
FRI-Ler allele accelerated flowering by 10 days in response 
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Figure 1 Functional comparison of the most common deleterious FRI alleles, (a) Leaf number at flowering from individuals in the Col x Ler 

recombinant inbred line set grouped by their genotype at FRI and FLC. (b) Vernalization response of Arabidopsis natural accessions measured as 
the reduction of days to flowering in plants vernalized for 5 weeks at 4°C compared to unvernalized plants. Accessions are grouped by their FRI 
alleles as described in [22]. The number at the bottom of each bar indicates the number of accessions in each group; error bars in (a) and (b) 
correspond to the standard error of the mean. Letters in each bar or point represent significance groups as determined by Tukey HSD test, (c) Exonic 
organization for FRI-wt (Bil-7), FRI-Col and FRI-Ler. Red regions indicate the indels characteristic for FRI-Col and FRI-Ler. (d) Predicted proteins for all three 
alleles shown in (c). The dashed region in FRI-Col represents a shift in the ORF caused by the 1 6 bp deletion. Shaded regions in the background 
indicate predicted coiled-coil domains. 



to vernalization (Figure lb). This acceleration was inter- 
mediate to the low response of accessions containing the 
FRI-Col deletion and the high response of accessions car- 
rying the FRI-wt allele (Figure lb). The same trend was 
found for the vernalization response quantified as total 
leaf number, although the difference between accessions 
carrying the FRI-Ler allele and the FRI-Col allele was not 
statistically significant (Additional file 1). 

Taken together, our results disagree with the common 
belief that the FRI allele present in Ler is not functional. 
However, from this analysis we cannot completely rule out 
the possibility that the observed phenotypes are caused by 
additional loci in linkage disequilibrium with FRL 

Structural characterization of the FRI-Ler and FRI-Col 
alleles 

We looked for further evidences of a functional FRI-Ler al- 
lele by studying sequence variation among Arabidopsis ac- 
cessions. We sequenced genomic DNA from FRI including 
its upstream and downstream regions in Col-0, Ler-1 and 
Bil-7, the latter being a winter-annual accession that con- 
tains a functional FRI-wt allele identical to the published 
H51 allele [2,14]. All polymorphisms found in FRI-Col and 
FRI-Ler correspond exactly to the ones described before, 
such as the insertion/deletions that define these two allelic 
classes (Figure Ic; [2,14]). The H51 allele is predicted to 
encode a fully functional 609 aa protein [GenBank: 



AAG23415] while the Col allele yields a 314 aa protein 
[GenBank: AEE81913], truncated due to a premature stop 
codon resulting from a 16 bp deletion in exon 1 (Figure Ic 
and d; [2]). In the case of FRI-Ler, a 376 bp deletion com- 
bined with a 31 bp insertion removes the translational start, 
but creates a new, out-of-frame start codon that is pre- 
dicted to yield a 41 aa protein [2]. Interestingly, in addition 
to the out-of-frame start codon, FRI-Ler contains an in- 
frame ATG codon downstream of the original start, which 
would result in a protein missing 42 aa of the N-terminus. 
Apart from this deletion and one conservative amino acid 
change from L to I in the first coiled-coil domain, the FRI- 
Ler protein is identical to the functional Bil-7 allele 
(Figure Id). There are no known functional motifs in the 
deleted segment of the FRI protein, and, thus, if transcrip- 
tion and translation occurred from the downstream in- 
frame start codon, the resulting FRI-Ler protein could pos- 
sibly be as functional as the full-length Bil-7 protein. 

We looked for evidences of transcription of this long 
FRI-Ler allele using RNA-seq reads from seven accessions 
containing the characteristic Ler indel [36]. We found 
signals of expression across the full length of the gene, 
including the C-terminal part, in all accessions (see 
Additional file 2). This raised the question of whether 
this long FRI-Ler transcript is translated into a protein 
and if so, whether it is sufficient to delay flowering and 
to confer vernalization sensitivity. 
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Comparing the effect of the FRI-Ler versus the FRI-Col 
allele in transgenics 

We studied the existence of a functional FRI-Ler allele by 
comparing its function with that of the FRI-Col allele in 
transgenic plants. For this, we cloned the putative coding 
region plus the upstream and downstream region of both 
alleles and transformed them in a common background. 
We chose Col-0 as a recipient because our previous results 
suggest that FRI-Col carries a loss-of-function mutation 
(Figure 1), and because Col-0 has been shown to contain a 
strong FLC allele that is activated in the presence of a FRI- 
wt allele [3]. First, we confirmed the lack of function of the 
FRI-Col allele by growing three independent T3 lines homo- 
zygous for the FRI-Col transgene {Col-0[FRI-Col]) under 
long day conditions in the greenhouse. In this experiment 
none of the transgenic lines was significantly different from 
wild type (Figure 2a). In a consecutive experiment in similar 
conditions we grew homozygous T3 transgenic lines con- 
taining the FRI-Ler transgene {Col-0[FRI-Ler]) and ob- 
served a significant delay in flowering in all lines compared 
with the Col-0 wild type (Figure 2b), confirming FRI-Ler 
functionality. 

Because FRI delays flowering time through upregulation 
of the floral repressor FLC [3], we performed qRT-PCR to 
test the expression of FLC in the transgenic plants. As ex- 
pected, we found elevated expression of FLC associated 
with the FRI-Ler transgene but not with the FRI-Col trans- 
gene (Figure 2c), despite all lines presenting increased FRI 
expression (Additional file 3). In addition, a prolonged 
cold treatment, known to abolish FLC expression, sig- 
nificantly reduced leaf number in Col-0[FRI-Ler] lines, 
but not in Col-0 wild type or in Col-0 [FT?/- Co/] lines 
(p < 0.001, p = 0.06 and p = 0.28 for the interaction 
between vernalization and genotype in a two-way ANOVA 
respectively, compare Figure 2d and Figure 2a and b). 

Altogether, our results demonstrate that the FRI-Ler 
allele, but not the FRI-Col allele, has the ability to delay 
flowering, and that this delay is achieved through upreg- 
ulation of FLC, 

Expression of FRI and FLC in natural accessions 

The functionality of the FRI-Ler allele is contrasting with 
its frequent presence among early flowering accessions 
[2,14]. Our analyses also show that, although accessions 
carrying FRI-Ler alleles have stronger vernalization re- 
sponses than accessions containing FRI-Col alleles, they 
do not reach the levels observed in lines carrying FRI-wt 
alleles (Figure lb and Additional file 1). This suggests a 
reduced functionality of the FRI-Ler allele in comparison 
with the FRI-wt allele, which can be caused by differ- 
ences in the protein sequence, in expression levels or in 
post-transcriptional/post-translational mechanisms. Pre- 
vious works have pointed to the low expression of the 
FRI-Ler allele, possibly as a consequence of the indel in 



0 

E 
c 

CD 



(a) 



30 



i 20 

o 



(b) 











a 


a 


a 


a 


8 


11 


11 


10 



(c) 



— I 1 1 1— 

0^1-1- 

I I I I 
o CO 10 

O 











rh 


c 




c 




11 


12 


12 



T 1 1 r 



O 1- ^ CM CO 

I I I I I 

TS ^ CD O 00 
CM 



C 

o 

CO ^ 
X ^ 

0 



20 
15 
10 
5 



Col-0 

Col-0 [FRI-Col\ 
Col-0 [FRI-Ler] 



a a a a 
3 3 3 3 



iimim 



be 



be 



(d) 



0 

C 0 

0 g 

H — • 

o 



20 - 
15 - 
10 - 
5 - 
0- 



11 



ab 



ab 



abc 



rh 

abe 



be 



o 



"I I I I I I r 



^ 1- 

-L I I 

O 1- CO 

O 



I 

10 



1- ^ CM CO 

I I I I 

^ CD O 00 

1- 1- CM 



Figure 2 Characterization of transgenic lines carrying the FRI-Col 
or FRI-Ler allele, (a) Flowering time expressed as total leaf number 
(rosette + cauline) from homozygous T3 single insertion lines containing 
a FRI-Col transgene are shown in comparison to the untransformed Col- 
0 wild type, (b) Experiment as in (a) but using homozygous T3 single 
insertion lines transformed with FRI-Ler. (c) Expression of RCfor the 
same genotypes as above. In this experiment, leaf tissue was collected 
from 10 day-old seedlings grown in long day conditions. The number in 
each bar indicates the number of biological replicates used. Expression 
was normalized to the expression of PP2A. (d) Flowering time after 
vernalization quantified as total leaf number for the same lines as 
above. Plants were vernalized for four weeks at 4°C and subsequently 
grown under long day conditions. The number in each bar indicates 
the number of individual plants per line analyzed. Error bars represent 
the standard error of the mean. Letters in each bar represent 
significance groups as determined by Tukey HSD test. 
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its promoter [2,37]. This led us to focus on a possible 
transcription inhibition as the cause for early flowering 
in natural accessions containing this allele. 

We analyzed the expression profiles of FRI and its 
downstream target FLC in 137 Arabidopsis accessions 
containing FRI-Ler, FRI-Col or FRI-wt alleles [38]. Acces- 
sions carrying FRI-Ler alleles showed significantly reduced 
expression of FRI when compared to accessions carrying 
FRI-Col or FRI-wt alleles, suggesting a regulatory defect in 
the FRI-Ler allele (Figure 3a). Interestingly, this residual 
FRI-Ler expression results in higher FLC levels than in the 
accessions carrying FRI-Col alleles, although this differ- 
ence is not significant (Tukey HSD test p value p = 0.079), 
possibly due to the large variation in FLC expression 
found among accessions with FRI-Ler alleles (compare 
error bars in Figure 3b). 

Allele-specific expression of FRI-Ler and FRI-Col 

There are two alternative causes for the lower FRI expres- 
sion found in Arabidopsis accessions carrying the FRI-Ler 
allele when compared to accessions carrying FRI-wt or FRI- 
Col alleles. Either the FRI-Ler allele contains a mutation in 
cis, such as the promoter indel that defines the class, or its 
reduced expression is caused by the effect of trans- 
regulators that may occur in the early flowering accessions. 
In order to distinguish among these possibilities, we ana- 
lyzed the allele specific expression of FRI using pyrose- 
quencing in hybrids from crosses of Arabidopsis accessions 
differentially expressing this gene. Because expression of 
both alleles in a hybrid cell is controlled by the same set of 
regulators, the relative abundance of each allele in the hy- 
brid should be equaled out if trans-regulators cause the dif- 
ferences between accessions. In contrast, if the expression 



differences between accessions are caused by mutations in 
cis, the differences in the expression of the alleles in the hy- 
brid will be maintained [39]. 

We generated Fi hybrids by crossing Ler-1 with Col-0 
and with four accessions containing FRI-wt alleles (Bor-4, 
Bur-0, Sha and Knox- 10; [14]). In all Fi individuals ana- 
lyzed the relative abundance of each allele measured in 
genomic DNA was close to 50%, which ensured that the 
assay had no preference for either of the two alleles 
(Figure 3c). Analysis of cDNA from the same individ- 
uals revealed lower relative abundance of FRI-Ler than 
either the FRI-wt or the FRI-Col alleles (Figure 3d). 
These differences suggest that at least part of the 
phenotypic differences found between accessions car- 
rying the FRI-Ler allele and those carrying the FRI-wt 
allele is due to a cis regulatory mutation reducing the 
expression of FRI-Ler. Interestingly, the differences in 
expression found among accessions classified by their FRI 
allele are larger than the differences observed in allele spe- 
cific expression in the hybrids (compare Figure 3a and d), 
suggesting the existence of additional trans-regulators 
controlling differences in expression among alleles. Fur- 
ther experiments will be required to determine the precise 
mode of regulation of this gene. 

Effect of coding polymorphisms in FRI-Ler 

In addition to the suggested regulatory polymorphisms 
causing inhibition of the expression of the FRI-Ler allele, 
its predicted protein differs from the FRI-wt protein in the 
deletion of 42 amino acids from the N terminal and in one 
substitution located in a coiled-coil domain (Figure Id). To 
test whether these coding differences have an effect on the 
fianction of the FRI-Ler protein compared to FRI-wt, we 
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Figure 3 Relative expression levels of FRI and FLC in Arabidopsis accessions and Fi hybrids grouped by their FRI allele, (a and b) 1 37 

Arabidopsis accessions were grown in long days under greenliouse conditions, sampled after four weeks and expression levels of FLC and FRI 
were quantified from northern hybridizations [38]. Letters in each bar indicate the significance groups as determined by a Tukey HSD test. The 
number in each bar indicates the number of accessions in each group, (c and d) Positions with SNPs specific to Ler-1 or Col-0 were targeted 
using pyrosequencing in genomic DNA (c) or in cDNA (d) of Fi hybrids. Leaf material for hybrids involving Bor-4, Bur-0, Knox-10 or Sha was 
collected at the time of bolting from plants growing under 12-hour photoperiods in an environmental chamber. Leaf material from hybrids 
involving Col-0 was collected from 10 day-old seedlings grown in the greenhouse in long days. Error bars indicate the standard error of the mean. 
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generated transgenic lines carrying constructs in which the 
coding region of each allele was placed downstream of the 
native promoter of the other allele (Figure 4). Due to con- 
straints (see Material and Methods), this experiment was 
performed using Ti lines carrying FRI-wt (the FRI-wt allele 
expressed under its native promoter), FRI-wt ::FRI'Ler (the 
FRI-Ler allele expressed under the FRI-wt promoter) or 
FRI-Ler::FRI-wt (the FRI-wt allele expressed under the 
FRI-Ler promoter) constructs, and individuals from two 
independent T3 lines carrying the FRI-Ler construct (the 
FRI-Ler allele expressed under its native promoter). As ex- 
pected, lines carrying the wild type promoter driving the 
expression of the wild type allele (FRI-wt) flowered with 
the same number of leaves as Col FRI-Sf2, a line in which 
the FRI-wt allele from the accession Sf-2 was introgressed 
into Col-0 [2]. In addition, lines carrying constructs with 
the FRI-Ler promoter flowered significantly earlier than 
those with alleles driven by the FRI-wt promoter. This 
supports our hypothesis that cis-regulatory elements spe- 
cific to the FRI-Ler promoter are the cause of early flower- 
ing in accessions carrying this allele. Interestingly, lines 
carrying the FRI-Ler coding region flowered earlier than 
lines carrying the FRI-wt coding region when expressed 
under the same promoter. Although this acceleration of 
flowering was not significant, it was observed both 
when using the FRI-Ler promoter and the FRI-wt pro- 
moter. This result leaves open the possibility that the 
loss of 42 amino acids and/or the mutation in the coiled 
coil domain in the FRI-Ler protein contribute to accel- 
erating flowering in the accessions that express it. 

Analysis of life history parameters in the FRI-Col and 
FRI-Ler accession groups 

A major question that derives from the classification of 
FRI-Ler as a semi-functional allele is whether it represents 
a change in the life history of the accessions that carry it 
when compared to the accessions carrying the FRI-Col al- 
lele. Although FRI-Ler accessions revealed a significant 
higher response to vernalization than accessions with the 
FRI-Col allele (Figure lb), both groups only include early- 
flowering accessions that are not likely to differ in their 
general life history strategy. Nevertheless, we investigated 
the variation in specific life history traits between the two 
allelic classes, which can result in a distinct ability to adapt 
to local environments. To address this question, we reana- 
lyzed fitness-related phenotypes for a panel of accessions 
grown at four different field locations [40]. Although we 
observed small fitness differences in accessions carrying 
the FRI-Col and FRI-Ler alleles at various locations, these 
were not significant (see Additional file 4). On the other 
hand, accessions containing the FRI-wt allele were often 
significantly different in fitness traits when compared 
with FRI-Ler and FRI-Col. Further experiments will be re- 
quired to determine if differences between the FRI-Ler 
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Figure 4 Flowering time of transgenic lines carrying 
promoter-swap constructs involving the FRI-Ler and FRI-wt 
alleles. Plants were grown in an environmental chamber set to long 
day conditions and the total number of leaves (rosette + cauline) 
were counted on the day of the opening of their first flower. Col-0 
is the untransformed control and Col FRI-Sf2 is a near isogenic line 
containing the FRI-wt allele from the accession Sf-2 [2]. Phenotypes 
were scored in Ti individuals for lines carrying the FRI-wt, FRI-wt::FRI-Ler 
and FRI-Ler.FRI-wt constructs, and in individuals from two independent 
T3 lines for lines carrying the FRI-Ler construct (see Material and 
Methods). Numbers in each bar indicate the number of individual 
plants analyzed per genotype. Error bars represent the standard error 
of the mean. Letters in each bar represent significance groups as 
determined by Tukey HSD test. 

V ) 

and FRI-Col allelic classes are relevant in specific environ- 
mental conditions (e.g. under stress). 

Discussion 

FRI, an activator of the negative regulator PLC, is the major 
determinant of natural variation in flowering time and 
vernalization sensitivity in Arabidopsis thaliana [3]. So far, 
natural allelic variants of FRI have been clustered into two 
groups, the functional FRI-wt group and the non-functional 
group that included the FRI-Col and the FRI-Ler alleles. 
Here we prove that FRI-Ler, the most common allele found 
so far in summer- annual accessions, is partially functional, 
although it is a weak allele compared to the FRI-wt allele 
present in winter-annual types. We suggest that the indel 
found in FRI-Ler, including part of the promoter and the 
beginning of the coding sequence, allows transcription of 
an active messenger from a downstream methionine. This 
resulting transcript shows lower expression and perhaps 
slightly reduced activity than the transcripts from both 
FRI-wt and FRI-Col, but has the ability to up regulate FLC, 
delay flowering time and confer vernalization sensitivity. 
In contrast, the FRI-Col allele is unable to confer any of 
these phenotypes, despite being expressed at a level similar 
to the FRI-wt alleles present in winter-annual accessions. 

Although FRI has a prominent role, not all variation 
in FLC expression and flowering time observed in Ara- 
bidopsis accessions can be explained by variation at this 
gene. Previous studies have demonstrated that allelic 
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variation at FLC itself provides a basis for the development 
of a summer-annual flowering time habit [24,25,41]. As 
shown by [15], variation in flowering time may be asso- 
ciated with FLC only in the presence of a functional FRI 
allele. Consistent with this, we observed that accessions 
carrying the FRI-Ler or FRI-wt alleles, but not those car- 
rying a FRI'Col allele, varied considerably in FLC ex- 
pression and flowering time (Figures 1 and 3). Further 
variation of FLC levels can be due to additional factors 
such as mutations in genes of the autonomous floral 
promotion pathway [25,42-44]. Moreover, allelic vari- 
ation at genes directly interacting with FRI, namely 
SUPPRESSOR OF FRIGIDA 4 (SUF4), FRIGIDA-LIKE 1 
(FRLl), FRIGIDA ESSENTIAL 1 (FESl) and FLC EXPRES- 
SOR (FLX), could have an impact [18]. For instance, it 
has been shown that Ler contains a non-functional 
FRLl allele, which is compensated by a functional FRI- 
GIDA-LIKE 2 (FRL2) allele [45]. The latter restores FRI- 
mediated up regulation of FLC expression and, thus, a 
late flowering phenotype. In contrast, Col-0 carries a 
functional FRLl allele able to interact with FRI, but an 
impaired FRL2 allele [45]. 

The effect on flowering time, FLC expression and 
vernalization response conferred by the FRI-Ler allele 
observed in the present study is not comparable to the 
effect of the FRI-wt allele. Although we show here that 
both alleles differ greatly in expression (Figure 3), we de- 
tect small, non-significant but consistent differences be- 
tween the FRI-Ler and the FRI-wt allele driven under the 
same promoter (Figure 4). This suggests that the puta- 
tive truncation of the first 42 aa of the N-terminal region 
and/or the amino acid substitution in FRI-Ler contribute 
to the functional difference between alleles (Figure 4). In 
fact, previous work demonstrated that a deletion of 118 
residuals from the FRI N-terminus resulted in reduced 
functionality. Nevertheless, plants missing that part of 
the protein still flowered later and revealed higher FLC 
expression levels than both Col-0 and transgenic lines 
lacking a substantial part of the FRI C-terminus [17]. 
Further studies used a deletion series to demonstrate 
that the ability of FRI to act as the scaffold protein for a 
FLC activating complex depends mainly on its C-terminal 
part [18]. These studies suggest that the truncated FRI-Ler 
protein should be able to interact with its known partners 
to build the FRI-complex and, thus, to induce FLC expres- 
sion. Finally, the amino acid change from leucine in FRI- 
wt to isoleucine in FRI-Ler in the first coiled-coil do- 
main is very conservative, but might impair the func- 
tionality of the protein (Figure Id). However, the small 
differences observed between the FRI-Ler and FRI-wt 
alleles expressed under the same promoter suggest that 
if real, the effect of these mutations on flowering time is 
minor compared to the effect of their differences in 
expression. 



Despite the functional difference between the two al- 
leles, accessions of the FRI-Col and FRI-Ler group are 
all rapid cycling and, thus, do not display a difference 
in their general life history strategy. Furthermore, dif- 
ferences in fitness traits between both allelic groups 
could not be detected in the present study. A more 
comprehensive analysis of larger data sets collected in 
multiple environments might be required to detect a 
distinct effect of the FRI-Ler allele in nature. In fact, 
[46] have demonstrated that, depending on the sea- 
sonal timing of specific environmental signals, Arabi- 
dopsis accessions in nature are capable of both life 
history strategies. For instance, Col-0 growing during 
autumn in the field at a specific location in Germany 
(Halle) displayed a winter-annual habit and flowered al- 
most at the same time as the introgression line Col 
FRI-Sf2 that carries a functional FRI-wt allele [46]. In 
contrast. Col FRI-Sf2 exhibited an unexpected rapid- 
cycling phenotype when grown in summer at one loca- 
tion in England (Norwich). Nevertheless, her resembled 
Col-0 with regard to flowering time at all locations 
tested [46]. Furthermore, we cannot exclude that the 
FRI-Ler allele has a role in the modulation of adaptive 
parameters other than flowering time and vernalization 
response. For example, pleiotropic effects of FRI have 
been reported for traits such as water use efficiency 
and leaf senescence, where the latter is closely linked to 
plant reproduction [47,48]. 

Conclusions 

In the present study we disprove the common assumption 
that the widespread FRI allele found in the Arabidopsis 
Ler accession is not functional. We demonstrate that the 
FRI-Ler allele increases flowering time and FLC expression 
and induces vernalization responses to levels that are in 
between the ones observed for the true null FRI-Col allele 
and the fully functional FRI-wt allele. 

These intermediate phenotypes observed in plants car- 
rying FRI-Ler could be explained by the presence of a 
close to full-length protein never before associated with 
this allele. We show that the differences in functionality 
between the FRI-Ler and the FRI-wt alleles are largely 
due to expression polymorphisms, although variation in 
the protein sequence may also play a role. 

Using the limited data available from plants grown under 
natural conditions, we cannot conclude that the FRI-Ler al- 
lele confers differences in fitness or life history strategies 
when compared to the FRI-Col allele. On the other hand, 
phenotypic differences between accessions carrying these 
alleles are clear and could be advantageous under specific 
natural environments not explored so far. 

In summary, we demonstrate the existence of an allelic 
series in FRL Besides completely functional and non- 
functional alleles, we have found a widespread allelic class 
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with intermediate functionality. The use of this novel clas- 
sification will increase the accuracy of adaptive studies in 
Arabidopsis. 

Methods 

Analysis of published datasets 

Leaf number from Recombinant Inbred Lines (RILs) in 
the Ler x Col population [31] was analyzed by grouping 
lines according to their genotypes at the closest molecular 
markers to FRI (m506, chromosome 4 at 0.0 cM) and FLC 
(g4560, chromosome 5 at 17.3 cM). Phenotypic and geno- 
typic data for these individuals was kindly provided by 
Johan W. van Ooijen and Caroline Dean. 

Responses to vernalization from 126 Arabidopsis acces- 
sions were obtained from published data [22] and analyzed 
by subtracting days to flowering with vernalization from 
days to flowering without vernalization (plants grown in 
controlled environment rooms with 16 hour light). In this 
dataset, FRI alleles are classified as FRI-Len FRI-Col, non- 
functional alleles or novel alleles. The rest of the alleles 
were assumed to be functional (FRI-wt) as they do not 
contain the Ler or Col-0 indels and the accessions that 
carry them are late flowering. For our analysis, we con- 
sidered only those accessions containing aUeles classi- 
fied as FRI-wt, FRI-Ler or FRI-CoL Accessions for which 
the FLC allele was classified as non-functional or novel 
were removed. 

Expression data for FRI and FLC was obtained from 
available repositories holding northern hybridization 
quantification for 137 Arabidopsis accessions (greenhouse 
grown, 16 h light [38]). Accessions in this dataset were 
classified according to their FRI alleles by combining ge- 
notypes from the following sources: [14,21,22,27,37]. The 
assignment of FRI alleles for each accession can be found 
in Additional file 5. 

Fitness data was obtained from Arabidopsis accessions 
grown in four different field locations [40]. These acces- 
sions were classified according to their FRI alleles as above. 

Sequencing of FRI in Col and Ler 

Genomic DNA of Col-0 and Ler-l was extracted from 
young leaves using a Plant DNeasy Mini Kit (Qiagen, 
Chatsworth, CA, USA). A region of approximately 3250 bp 
including the complete FRI gene (AT4G00650) and an up- 
stream and downstream region was amplified in overlap- 
ping fragments of 600-700 bp using Phusion High-Fidelity 
DNA Polymerase (New England Biolabs). As a reference, 
the same region was sequenced for the accession Bil-7 car- 
rying a fully functional FRI-wt allele [14]. Pooled PCR prod- 
ucts from four independent reactions were purified using 
QIAquick PCR purification kit (Qiagen, Chatsworth, CA, 
USA) and sequenced via Sanger sequencing at the Max 
Planck Genome Centre Cologne. The individual sequences 



were assembled and aligned against the sequence of Bil-7 
using SeqMan Pro (DNAstar, Madison, WI, USA). 

Cloning and phenotyping of transgenics 

In order to clone the FRI-Ler, FRI-wt and FRI- Col alleles, 
we designed primers flanking positions -1372 to +2691 
relative to the FRI start codon annotated in TAIRIO, 
which include the complete gene plus an upstream region 
of 1061 bp and a downstream region of 379 bp. This re- 
gion was amplified by PCR from DNA of the accessions 
Ler-l, Bil-7 and Col-0 using Phusion High-Fidelity DNA 
Polymerase (New England Biolabs). For FRI-Ler and 
FRI-Col, PCR fragments were introduced into the binary 
vector pCAMBIA2300 making use of EcoRI/BamHI re- 
striction sites. For FRI-wt, the PCR fragment was intro- 
duced into binary vector pBinGlyRedl (kindly provided 
by Ed Cahoon, University of Nebraska) making use of 
the same EcoRI/BamHI restriction sites. 

Promoter-swap constructs were generated using the 
MultiSite Gateway® Pro 2.0 system (Life Technologies) 
and genomic DNA from the accessions Ler-l {FRI-Ler) 
and Bil-7 {FRI-wt), The FRI-Ler promoter region was 
cloned starting from position -1372, relative to the an- 
notated FRI start codon in TAIRIO, and ending at pos- 
ition + 126, which is right before the first ATG 
downstream of the indel described in FRI-Ler, For FRI- 
wt, the promoter region cloned ranged between posi- 
tions -1372 and -1. Both promoter sequences were in- 
troduced into pDONRTM 221 Pl-P5r by BP reaction. 
Simflarly, the coding sequences for both alleles to- 
gether with 379 bp of its downstream sequence (pos- 
ition + 1 to position +2691 bp for FRI-wt and position 
+ 127 to position +2691 for FRI-Ler) were cloned and 
introduced into pDONRTM 221 P5-P2. Finally, FRI- 
wt::FRI-Ler and FRI-Ler::FRI-wt constructs were gener- 
ated from LR reactions using the binary vector pFAST. 

All final constructs were transformed into E, coli strain 
ToplO (One Shot® TOPIO chemically competent E. coli, 
Invitrogen). Inserts ft-om positive colonies were sequenced 
and verifled by comparing them to the PCR template se- 
quence or to the expected in silico constructs. Subsequently, 
aU constructs were transformed into Agrobacterium 
tumefaciens strain GV3101 and used for transformation 
of Col-0 plants by the floral dip method [49] . 

Transgenic plants carrying FRI-Ler or FRI- Col constructs 
(in pCAMBIA) were selected on MS agar plates containing 
50 mg/1 kanamycin. Independent Ti plants carrying a single 
insertion were identified and homozygous T3 plants were 
used for subsequent analyses. Flowering time in T3 lines 
carrying FRI-Col or FRI-Ler was measured in the green- 
house under long day conditions (16 h day length) with 
and without vernalization treatment for four weeks at 4°C. 

Ti plants carrying the promoter swap between FRI-Ler 
and FRI-wt (in pFAST) and those carrying the Bil-7 



Schmalenbach et al. BMC Plant Biology 2014, 14:218 
http://www.bionnedcentral.conn/1471-2229/14/218 



Page 9 of 1 1 



FRI-wt alleles (in pBinGlyRedl) were selected by their 
red fluorescence under a fluorescence stereomicroscope 
(Leica MZ16F, Leica Microsystems, Germany) using 
green light of wavelength about 580 nm. Flowering time 
was recorded in Ti lines carrying the FRI-wt allele and 
the promoter swap lines FRI-wt:FRI-Ler and FRI-Ler:: 
FRI-wt under long day condition in an environmental 
chamber. This experiment included plants from two in- 
dependent T3 lines carrying the FRI-Ler::FRI-Ler allele 
(lines 14-1 and 16-4 in Figure 2). Using these T3 lines 
allowed us to sow all the plants in the experiment dir- 
ectly on soil, as FRI-Ler Ti lines do not have red fluores- 
cent marker and would have required selection on 
kanamycin plates. 

In all experiments, flowering time was quantified as 
total leaf number (rosette + cauline) on the day the first 
flower opened. 

Expression analysis using quantitative RT-PCR 

For expression analysis of transgenics, seedlings (leaves plus 
shoot) grown in an environmental chamber at long day 
conditions (16 h light) were harvested 10 days after sowing. 
Material of 8 to 10 seedlings was pooled for each of three 
biological replicates. RNA was extracted using Trizol 
(Ambion® TRIzol® RNA Isolation Reagent, Life Technolo- 
gies) and transcribed into cDNA using Super Script® II Re- 
verse Transcriptase (Invitrogen). Quantitative RT-PCR was 
performed on a CFX384 Touch™ Real-Time PGR Detection 
System using SYBR Green dye (iQ™ SYBR® Green Super- 
mix, Biorad) using the following gene-specific primers: 
FIC-fwd: 5 -CCGAACTCATGTTGAAGCTTGTTGAG- 
3', FIC-rev: 5 - CGGAGATTTGTCCAGCAGGTG-3', 
m-fWd: 5 - TGCCTGATCGTGGTAAAGGGAAG-3' and 
FRI-vew: 5 - AGCACCGGCAATCTCATTCGAAC-3'. Ex- 
pression values were determined using the standard curve 
method and normalized to the expression of PP2A {PP2A- 
fWd: 5'-TAACGTGGCCAAAATGATGC-3', PP2A-rev: 5 - 
GTTCTCCACAACCGCTTGGT-3'). Normalized expres- 
sion was averaged for three biological replicates each ana- 
lyzed in two or three technical replicates. 

Allele-specific expression analysis 

Leaves from Fi hybrids involving the Ler-1 parent were col- 
lected at bolting from plants grown in 12 h days in a con- 
trolled environmental chamber. Material from F^ hybrids 
involving the Col-0 parent were collected from ten-day old 
seedlings grown in long days (16 h light) in a controlled en- 
vironmental chamber. Here, material of 8 to 10 seedlings 
was pooled per biological replicate. RNA was extracted 
from 20 mg of frozen material using RNeasy Plant kit 
(Qiagen) combined with on-column DNase digestion using 
RNase-Free DNase Set (Qiagen). Subsequently, cDNA was 
synthesized using Superscript® III Reverse Transcriptase 
kit (Invitrogen) in the presence of RNase inhibitor 



RNasin (Promega) following specifications from the 
manufacturer. Genomic DNA was extracted from the 
same samples using DNease Plant kit (Qiagen). 

PGR distinguishing FRI-Ler versus FRI-wt alleles tar- 
geted a variant from C (in Bor-4, Bur-0, Sha and Knox- 
10) to A (in Ler-1) at chromosome 4, position 269257 
(TAIRIO). The fragment containing this SNP was ampli- 
fied using primers 5'Biotin-TCAGTTGCAGTGGAAA 
CATTCA-3' and 5'-GCGTTTTCGATTGACTCGATG 
T-3', and pyrosequencing was performed using primer 
5'-TGACTCGATGTGCTTCT-3'. To distinguish the 
FRI-Ler and FRI-Col alleles in the Ler x Col F^ hybrid 
we targeted a variant from A in Col-0 to G in Ler-1 at 
chromosome 4, position 269469 (TAIRIO). The frag- 
ment containing this SNP was amplified using primers 
5 'Biotin-ATTGTACCGGAGACGTCGAATAA-3 ' and 
5'-GGCCAATTTCAAAGCTGAAG-3', and pyrose- 
quencing was performed with primer 5'-CTTTGCTA 
CACATCAACTC-3'. 0.5 \iL of cDNA was used for 
PGR, which was conducted in a solution (2.5 units of 
Taq polymerase, 2.5 mM MgCl2, 200 \iM dNTPs, buf- 
fer B, 0.2 mM of each dNTP; Bio-Budget) containing 
20 pmol of each primer, in the total volume of 25 [xL, 
The PGR products were obtained with 50 cycles (93°G, 
45 sec; 60°G, 45 sec; 72°G, 1 min). 

Allele-specific mRNA abundances were measured using 
PyrosequencerAB (Biotage AB, now Qiagen) following 
manufacturers instructions for sample preparation and 
pyrosequencing reactions. Vacuum sample preparation 
was performed using 15 (iL of PGR product mixed with 
5 [iL of Streptavidin Sepharose beads (GE Healthcare), 
40 [iL of PyroMark Binding Buffer (Qiagen) and 20 [iL 
of LiGhroSolv water (Merck). Pyrosequencing was per- 
formed in a PyroMark Q96 Plate Low (Qiagen) in each 
well containing 1 [iL of sequencing primer (10 [iM) and 
40 [iL of Pyromark Annealing Buffer (Qiagen) using 
PyroMark Gold Q96 Reagents (Qiagen). Allele-specific 
expression for each SNP was estimated using the pyrose- 
quencing software (PSQ 96MA 2.1.1, Biotage AB) based 
on the peak height for each allele at the SNP. A peak cor- 
rection factor of 0.86 was used for incorporation of dATP 
aS, as recommended by the manufacturer. 

Availability of supporting data 

The coding DNA sequence and translated protein sequence 
of the FRI-Ler allele supporting the results of this article are 
available through NGBIs GenBank under accession num- 
ber KJ545576 (http://www.ncbi.nlm.nih.gov/genbank). 

Additional files 



Additional file 1: Response to vernalization quantified as the 
decrease in the total number of leaves at flowering for Arabidopsis 
accessions grouped by their FRI allele, as described in [22]. 
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Additional file 2: RNA-seq coverage at the FRI locus in seedlings of 
seven Arabidopsis accessions carrying the FRI-Ler allele. 

Additional file 3: Expression of FRI in transgenic lines carrying the 
FRI transgene. 

Additional file 4: Assessment of fitness differences between 
accessions containing FRI-Col, FRI-Ler or FRI-wt alleles. 

Additional file 5: Assignment of FRI alleles in Arabidopsis 
accessions used in Figure 3 and Additional file 4. For both articles 
referenced in the table [38,46], FRI alleles are marked with y (yes) when 
identified in the article and with n (no) when not present in the article. 
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